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An HST parallax of the distant cataclysmic variable V1223Sgr, 
its system parameters, and accretion rate* 
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Abstract. Using the Hubble space Telescope Fme Guidance Sensor, we have measured the trigonometric parallax 
of the bright cataclysmic variable 1223 Sgr. The absolute parallax is TTabs = 1.96 ± 0.18 mas, making V1223 Sgr 
the most distant CV with a well-determined trigonometric parallax. This distance, a Lutz-Kelker correction, 
and the previously measured extinction yield an absolute visual high-state magnitude My = 4.0 ± 0.2. We 
outline a model, which is consistent with the observed spin-down of the white dwarf and provides for much of 
the UV/optical emission by reverberation of X-rays. From previous X-ray and UV/optical data, we derive an 
accretion luminosity iacc ~ (2.6 ± 0.8) 10^* erg s^^ , a white dwarf mass A4i = 0.93 ± 0.12 Mq, and an accretion 
rate M = (1.4 ± 0.3) 10^^gs~\ 
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1. Introduction 

V1223Sgr (4U1849-31) is one of the brightest confirmed 
intermediate polars, both optically and in X-rays. The 
only previous distance estimate of 600 pc was based 
on the observed E'b-v = 0.15 and a mean reddening of 
0.25 mag kpc"^ (Bonnet-Bidaud et al. 1982). Since V1223 
Sgr is not included in ground-based parallax programs, we 
have obtained a high-precision parallax using the Hubble 
Space Telescope Fine Guidance Sensor (FGS). This obser- 
vation increases the number of CV parallaxes measured 
with the HST FGS to eight (Harrison et al. 1999, 2000, 
2003, McArthur et al. 1999, 2001, Beuermann et al. 2003). 

V1223Sgr shows modulations at the orbital period 
^oib=3.366hrs, the spin period of the white dwarf 
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-Pspin=745.5s (Osborne et al. 1985), and the beat period 
Pboat=794.38s (Steiner et al. 1981). The UV/optical spec- 
tral energy distribution was interpreted by Bonnet-Bidaud 
et al. (1982) and Mouchet (1983) in terms of an accretion 
disk which reaches close to the white dwarf and should 
cause it to be spun up. Unexpectedly, however, the white 
dwarf was found to spin down (Jablonski & Steiner 1987, 
van Amerongen et al. 1987) suggesting that the disk stays 
much further away from the white dwarf. Our accurate dis- 
tance allows us to determine the luminosity of V1223 Sgr, 
estimate its accretion rate, and to show that a standard 
intermediate-polar model is consistent with the available 
observations. 



2. Observations, data reduction, and results 

2.1. HST observations 

The process for deriving a parallax for a cataclysmic vari- 
able from HST FGS observations has been described in 
papers by McArthur et al. (2001, 1999) and Harrison et 
al. (1999). The process used here is nearly identical to 
those efforts, as well as our recent program on EX Hya 
(Beuermann et al. 2003). An FGS program consists of a 
series of observations of the target of interest, and a set of 
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Table 1. Photometric and spectroscopic data for the V1223 Sgr reference frame. 



Star 


V 


B-V 


V-R 


V-I 


J-H 


H-K 


K 


Ay (mag) 


Spectral Type 


V1223 Sgr 


13.18 


0.02 


0.05 


0.15 


0.16 


0.06 


12.68 


0.47^ 




Ref-#1 


12.81 


1.80 


1.04 


2.24 


1.08 


0.14 


8.07 


0.81 


MIIII 


Ref-#2 


12.34 


0.59 


0.33 


0.66 


0.24 


0.11 


10.89 


0.41 


F5V 


Rcf-#3 


13.01 


1.16 


0.59 


1.13 


0.57 


0.12 


10.40 


0.40 


G9.5III 


Ref-#4 


13.04 


1.17 


0.60 


1.18 


0.70 


0.09 


10.28 


0.50 


G9.5III 


Ref-#5 


13.44 


1.40 


0.71 


1.34 


0.70 


0.15 


10.30 


0.84 


Kim 



^ Prom £b-v = 0.15 (Bonnet-Bidaud et al. 1982). 



Table 2. Astromcrtric; data and dca'iv(xi spcK;troscopic parallaxc^s for th(^ V1223 Sgr rcfcr(ui('(^ fram(\ 



Star 


02000 


PMa (mas/yr) 


^2000 


Y>Ms (mas/yr) 




n 






TT 






UCAC2 


FGS 




UCAC2 


FGS 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


V1223 Sgr 


18 55 02.3 


+2.2±2.6 


+2.0±0.2 


-31 09 49.5 


-23.7±2.9 


-24.6±0.2 


-67.090 


790.152 


0.37 


0.35 




Ref-#1 


18 54 53.3 


-8.1±3.8 


-9.2±0.2 


-31 10 21.4 


- 8.1±4.4 


- 5.8±0.2 


50.571 


809.254 


0.45 


0.36 


0.33 


Ref-#2 


18 54 42.3 


-0.2±2.6 


-0.5±0.2 


-31 09 10.1 


-13.7±3.0 


-13.9±0.2 


183.174 


722.876 


0.60 


0.49 


1.91 


Ref-#3 


18 54 49.2 


-7.9±3.8 


-5.5±0.2 


-31 07 59.6 


-18.0±4.4 


-17.8±0.2 


87.710 


662.492 


0.42 


0.45 


0.40 


Ref-#4 


18 55 01.5 


+1.3±3.8 


-0.7±0.2 


-31 07 53.2 


- 7.4±4.4 


- 7.4±0.2 


-70.136 


673.396 


0.41 


0.45 


0.41 


Ref-#5 


18 55 07.1 


+3.3±3.8 


+4.4±0.2 


-31 09 55.7 


- 2.4±4.4 


- 0.9±0.2 


-128.303 


803.079 


0.49 


0.54 


0.37 



^ ^ and r? are the relative positions in arcseconds from Eqs. (3) and (4), and a,, their standard deviations. 



four or more reference stars located close to that target. 
Typically, three epochs of observations, each comprised of 
two or more individual if5T pointings (orbits), are used 
to solve for the variables in the series of equations that 
define a parallax solution. For V1223 Sgr, wc obtained 
observations on four different epochs (2000 September, 
2001 March, 2001 September, and 2002 September) during 
which nine orbits of HST time were used. The extra epoch 
of observation was essential in deriving the high precision 
in the smallest parallax we have yet measured using the 
FGS. Extensive calibration data, as well as estimates of 
the distances and proper motions of the reference stars, 
are required to obtain a robust parallax solution. 

2.2. Spectroscopic parallaxes of the reference frame 

We have used a combination of spectroscopy and pho- 
tometry to estimate spectroscopic parallaxes for the ref- 
erence stars. The optical BVRI photometry of the five 
reference stars was obtained on 2001 March 13 using the 
CTIO 0.9 m telescope and the Cassegrain Focus CCD im- 
ager^. These data were calibrated to the standard system 
using observations of Landolt standards. The final photo- 
metric data set is listed in Table 1. Included is the Two 
Micron All-Sky Survey (2MASS) JHK photometry of the 
reference stars, transformed to the homogenized system of 
Bessell & Brett (1988) using the transformation equations 
from Carpenter (2001). Typical error bars on the photom- 
etry are ± 0.02 mag for the y-band measurements, and ± 
0.03 mag for the optical colors. The 2MASS photometry 
for the five reference stars has error bars of ± 0.02 mag. 



Optical spectroscopy of the reference stars, and a num- 
ber of MK spectral type templates, was obtained on 2001 
March 9 and 10 using the CTIO 1.5 m telescope with the 
Cassegrain Spectrograph^ with the 831 //mm "G-47" grat- 
ing and a two arcsecond slit, giving a spectral resolution of 
0.56 A/pix. We estimated spectral types of the reference 
stars from a comparison of the spectra with those of the 
MK-templates and present the results in the last column 
of Table 1. The visual extinctions of the reference stars 
(column 9 of Table 1) were estimated from a comparison 
of their spectral types and the photometric data, using the 
standard relations from Reike & Lebofsky (1985). Finally, 
spectroscopic parallaxes were derived using the spectral 
types and visual extinctions. The results arc listed in the 
last column of Table 2. To determine these values, we used 
the Hipparcos calibration of the absolute magnitude for 
main sequence stars by Houk et al. (1997), and that for 
giant stars tabulated by Drilling & Landolt (2000). The 
distinction between main sequence stars and giants is fa- 
cilitated by (i) the dichotomy of the My-distribution in 
the Hertzsprung- Russell diagram for stars of spectral type 
G and later and (ii) the ability of the Bayesian-like astro- 
metric software to identify an object with a true parallax 
substantially different from the one used as the starting 
value in the iterative approach. Details of the procedure 
are given in our previous papers cited above. For the as- 
trometric solution discussed below, we assumed error bars 
of ± 20% on our spectroscopic parallaxes. 

The reference stars span a wide range in spectral type, 
and have significant extinctions. Four of the five reference 
stars turn out to be late- type giants, providing the excep- 



^ http://www.ctio.noao.edu/cfccd/cfccd.html 



^ http://www.ctio.noao.edu/spectrographs/60spec/manual/ 
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tionally quiet reference frame, which is an essential for the 
high quality of the parallax derived for V1223 Sgr. Three 
of them (#2, #3, and #4) have extinctions similar to that 
of V1223Sgr which has Ay = 0.47 ±0.02 (Bonnet-Bidaud 
et al. 1982). It is somewhat disquieting that Ay does not 
correlate better with tt. A similar situation, however, was 
also found in the fields of RRLyr (Benedict et al. 2002a) 
and 6Cep (Benedict et al. 2002b). 

2.3. The astrometric solution 

The data reduction process for deriving a parallax of the 
target from FGS observations was identical to previous 
efforts, except for the fact that the astrometer has changed 
from FGS3 to FGSIR, successfully calibrated by Benedict 
et al. (2001). With the positions measured by FGSIR, we 
determine the scale, rotation, and offset plate constants 
for each observation set relative to an arbitrarily adopted 
constraint plate. The solved equations are 



y = 
V = 



x + lc[B-V) (1) 

y + lc{B - V) (2) 

Ax' + By' + C + R^{x'^ + y'^) - - P„^, (3) 

Dx' + Ey' + F + Ry{x'^ + y") - Mytrci - (4) 



where x and y are the rectangular HST coordinates and 
Ic is the derived lateral color correction for color B — V 
(Benedict et al. 1999). A. B. D and E are a set of scale 
plate constants, C and F are zero point corrections, 
and Ry are radial terms, and /ij, arc the proper motions, 
iroi is the time of the observation minus the constraint 
plate time, Pa and Pg are parallax factors, and and tTj, 
are the parallaxes in x and y. 

As in all our previous astrometric analyses, we em- 
ployed GaussFit (Jefferys, Fitzpatrick, & McArthur, 1987) 
to minimize . Since we had five high-quality reference 
star observations at each epoch, we were able to use 
the eight-parameter model, instead of the six-parameter 
model used in our previous parallax measurements, which 
had A ^ E and B = —D. The additional parameters 
yielded a significantly improved compared with the 
simpler model. The input proper motions for V1223Sgr 
and the reference stars were taken from the USNO 
CCD Astrograph Catalog, the "UCAC2" (Zacharias et 
al. 2004)^^, and are listed in columns 3 and 6 of Table 2. 
The Bayesian-like reduction software models the proper 
motions and the spectroscopically determined parallaxes 
of the reference stars as observations with errors to pro- 
duce an absolute (not relative) parallax for V1223Sgr. It 
yields the final proper motions listed in columns 4 and 7 
of Table 2. The internal consistency of the solution ob- 
tained for the V1223 Sgr field is very good. Fig. ^ shows 
the histograms of the residuals for all 220 individual x 
and y positional measurements of the target and the refer- 
ence stars. Gaussian fits give standard deviations of 0.88 
and 0.90 mas, respectively, which is better than usually 



obtained for parallax measurements with the HST FGS. 
Columns 8 to 11 of Table 2 list the "catalog" positions 
^ and rj and their standard deviations ct^ and Un which 
result from the fit of Eqs. (3) and (4) to the time series of 
the individual observations. The frequency- weighted mean 
values of these standard deviations are — 0.43 mas and 
cr,, ~ 0.46 mas. Reference star #3 has standard deviations 
in both coordinates as low as the other stars, confirming 
that its input parallax is not grossly in error. Hence, we 
must presently accept that it is a giant at a distance of 
^ 2.5 kpc, although its implied tangential velocity is large. 

The final absolute parallax of V1223 Sgr is tt = 1.96 
± 0.18 mas. The fine precision of this result is directly 
attributed to the very stable reference frame afforded by 
the distant giants. The limiting factor in the error budget 
turns out to be the remaining uncertainty in the proper 
motions of the reference frame. 

2.4. The distance to V1223 Sgr 

The measured distance to V1223Sgr is d ~ 1/tt = 
5IOI43 pc. As pointed out by Lutz & Kelker (1973), paral- 
laxes suffer from a systematic error in addition to the ob- 
servational one in the sense that the most probable true 
parallax is smaller than the observed absolute parallax 
because the number of objects per 7r-interval increases as 
TT^" with n = 4 for a constant space density. More de- 
tailed analyses of the problem including the magnitude 
and space density distributions of the parent population 
have been considered, e.g., by Hanson (1979), by Smith 
(1987), and by Oudmaijer et al. (1998). Here, we have sim- 
ply assumed that the parent population has an absolute 
magnitude at the brighter end of the range < My < 10, 
giving an index of n = 3.0 for the space density (Hanson 
1979). The Lutz-Kelker correction is usually quoted as a 
correction in visual magnitude and, for this choice of pa- 
rameters, is AMy — —0.071. The most probable distance 
modulus of V1223 Sgr including the Lutz-Kelker correc- 



tion is m — il/ = 8.61 



+0.21 

-0.19 



and the most probable dis- 



ftp://cdarc.u-strasbg.fr/cgi-bin/VizieR, catalogue 1/289 



tance including this correction is d = 527_43 pc. This is 
the largest distance of a CV so far measured with the HST 
FGS. 



3. System parameters and accretion rate of 
V1223 Sgr 

3.1. Absolute magnitude 

In its normal high state, V1223 Sgr has a long-term mean 
magnitude oi B — 13.1 (Garnavich & Szkody 1988), 
around which it varies by ±0.5 mag. Of this, the orbital 
and spin modulation accounts for about ±0.3 mag. There 
are also long-term trends and a time span between 1937 
and 1951 when it experienced drops to much fainter states 
(Garnavich & Szkody 1988). With B - y ~ in the 
high state (Table 2), V and B are practically the same. 
From the depth of the 2200A feature, Bonnet-Bidaud et 
al. (1982) derived E'b-v = 0.15 with an error of about 
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Table 3. Integrated observed and inferred fluxes for indi- 
vidual wavelength ranges (see text). 



Fig. 1. Histograms of the x and y residuals obtained from 
modeling V1223 Sgr and its reference frame. 

0.01, which yields a visual absorption Ay = 0.47 ± 0.03. 
This value of the absorption, the mean visual magnitude 
of = 13.1, and our new distance modulus give an ab- 
solute mean visual magnitude in the normal high state of 
My = 4.0 ± 0.2. This value has still to be corrected for 
inclination effects. 

3.2. Secondary star and inclination 

Penning (1985) and Watts et al. (1985) assumed that the 
mass of the secondary in V1223Sgr is that of a Roche- 
lobe filling main sequence star, M2 = 0.40 Mq. In many 
CVs, however, the secondaries seem to be expanded over 
their main sequence radii and have smaller masses than 
the main-sequence assumption suggests (e.g. Beuermann 
et al. 1998). We adopt M2 = 0.40 as the nominal 
secondary mass and comment on the effect a mass as low 
as 0.25 Mq will have on the the results. For a spectral 
class M4, a secondary in this mass range has a ii'-band 
magnitude of = 15.1 — 15.4 and contributes 8-11% of 
the observed JsT-band flux (Table 1). 

According to Watts et al. (1985), the radial velocity 
amplitude of the orbital motion of the white dwarf as 
derived from the wings of the Balmer emission lines is 
Ki = 56 ± 13kms~^. Penning (1985), using a superior 
time resolution, found that the Balmer and HeIIA4686 ra- 
dial velocities were modulated, in addition, at the spin 
period (see his Figs. 9 and 3 and note that the spin modu- 
lation is clearly present in his Fig. 3) . Pcnning's Ki agrees 
closely with that of Watts et al., which we adopt here. 
As noted already by these authors, the implied inclina- 
tion is small. For Mi = 0.93 Mq (see below), wc obtain 
i = (24 ± 7)° with M2 = 0.40 M0 and i = (37± 10)° with 
M2 = 0.25 Mq. 



3.3. X-ray and XUV emission of V1223Sgr 
3.3.1. The X-ray spectral analysis of Beardmore et al. 



Band 



Component 



Flux 
(erg cm~^s~^) 



0.0136-100 keV 

> 54 eV 
13.6 -54 eV 
912 - 1250 A 

> 1250 A 



total emitted X-ray flux^ 
additional XUV source 
not observable 
unobserved FUV^ 
observed UV/optical^/IR'* 



3.30 10"^° 
< 10-" 

? 

(2 - 3) 10-^° 
7.59 10"^° 



^ Integral of 'incident' spectrum of Beardmore et al. (2000); ^ 
extrapolated from A > 1250; ^ total mean flux from Mouchet 
(1983); " this work. 

analysis of Beardmore et al. (2000). These authors fit- 
ted a multi-temperature shock model to the ASCA and 
Ginga X-ray spectra of V1223 Sgr, which included a treat- 
ment of the X-albedo and of internal absorption within 
the source. From this fit, they deduced a shock tem- 
perature of kTg = 43I12 keV and the emitted spectrum 
fcmi{E) as a function of photon energy E (called 'in- 
cident' spectrum by them). Integration of fe-m\{E) over 
all photon energies above the Lyman edge gives the to- 
tal fiux = 3.30 10-1° erg cm-^s-i (Table 3, line 1), 
which we call 'X-ray flux' because only 3% of it is emit- 
ted below 0.1 keV. Fx refers to kTs = 43keV and scales, 
to a first approximation, as kTg. The X-ray luminosity Lx 
then scales as kTg, too. The emitted spectrum fg^i[E) 
is based on the observed ASCA spectrum for photon en- 
ergies E > 0.4 keV. but for still lower energies it merely 
represents an extrapolation of the multi-temperature ther- 
mal model spectrum. Wc demonstrate in the next Section, 
however, that a substantial additional XUV source does 
not exist in V1223 Sgr. 

The observed X-ray luminosity based on the 
Beardmore et al. (2000) emitted spectrum and our 
Lutz-Kelker corrected distance from Sect. 2.4, is Lx,i= 
4TT(fF^ = (1.11 ±0.36) lO^^ergs-S where the error ac- 
counts for the uncertainties in the normalization of the 
observed X-ray flux, in kTg, and in d. We assume that Lx,i 
represents the luminosity of one accreting pole and that 
the second pole is hidden behind the white dwarf. Such 
a geometry is generally expected if the shock height h is 
small compared with the white dwarf radius i?i , as it is in 
the model of Beardmore ct al., and still holds for a finite 
shock height if the inclination is small and the magnetic 
axis is more or less aligned with the rotational axis. We 
account for the second pole by writing Lx = (1 + /?) -^x,i 
and use /3 = 1, below. 



3.3.2. The quest for additional XUV emission 



In magnetic CVs of the polar subtype, a large fraction 
of the accretion energy is released as soft X-rays and 
Our derivation of the white dwarf mass Mi and the ac- dominates the total X-ray flux, which is used to estimate 
cretion rate M rests on the results of the X-ray spectral their accretion rates (Beuermann & Burwitz 1995). Does a 
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soft X-ray component exist also in the intermediate polar 
V1223Sgr? 

The XUV luminosity can be estimated by the Zanstra 
method using the observed line flux ^4586 of HcIIA4686 
photons with hv = 2.65 eV in ergcm~^s~^. The number 
flux of ionizing photons needed to produce this line is 



N- 



>54 



^4686 



Q!B 47r 
£4686 ^ 



(5) 



where U is the solid angle subtended by the photoionizcd 
gas, aB is the total recombination coefficient of Hell for 
the Baker & Menzel Case B in cm-^s"^, and £4686 the emis- 
sivity in the A4686 line in ergcm^s"^. For the pure recom- 
bination case, £4686/111^4686 is the effective recombination 
coefficient for this transition. We use the results of Storey 
& Hummer (f 995)^, which take full account of coUisional 
processes and are tabulated for electron number densities 
of Ne = 10^ — 10^^ cm~^ and a wide range of electron tem- 
peratures. They yield £4686/(h'^4686 ae) = 0.216 (O.f 40), 
0.157(0.125), and 0.210(0.156) for = 10^, lO^o, and 
10^^ cm~^ and electron temperatures Te = 10^ (5 10^) K, 
respectively. The particle density in the photoionizcd re- 
gion of V1223 Sgr may range from > lO^** cm"'^ in the 
disk and in the densest parts of the magnetically con- 
fined accretion stream to much lower values elsewhere in 
the magnetosphere or in the corona of the disk and its 
bulge, while the temperature is likely between 10^ and 
a few times 10^ K. The above numbers demonstrate that 
£4686 / (h'^4686 Q^B ) = 0.18 ± 0.05 is appropriate for a wide 
range of parameters in the photoionizcd regions in V1223 
Sgr. 

Steiner et al. (1981) quote an energy flux of the A4686 
line in V1223Sgr of 9 10~^^ ergcm~^s~^ and other pa- 
pers suggest values within a factor of two in both direc- 
tions (Bonnet-Bidaud et al. 1982, Penning 1985, Watts 
et al. 1985). We adopt the Steiner et al. value which 
becomes -^4686 = 1-6 10""'^^ crgcm~^s~-'^ after correction 
for reddening. Because of the large internal absorption at 
soft X-ray energies, we use fl ~ 47r. This choice of pa- 
rameters yields A^>54 ~ 0.21 photons cm~^s~'^, with an 
uncertainty of about a factor of two from the error in 
F4686 and the uncertainties in O and alfgg/aB- For com- 
parison, the Bcardmorc ct al. X-ray spectrum /cmi(-E') 
provides 0.11 photons cm~^s~^ with energies above 54eV, 
of which 50% are below 0.3 keV. As argued above, this 
number probably refers only to the pole facing the ob- 
server. The number from both poles would then be 
0.22 photons cm~^s~^, in agreement with the number 
needed to produce the observed A4686 flux. We conclude 
that there is no evidence for an additional component of 
XUV photons. The conservative limit to the energy flux 
of a potential XUV source given in line 2 of Table 3 refers 
to a source which provides the same number of ionizing 
photons as the X-ray source. 



ftp://cdarc.u-strasbg.fr/cgi-bin/VizieR, catalogue VI/64 



3.4. An intermediate polar model for V1223Sgr 

For further discussion of V1223Sgr, we adopt a simple 
intermediate polar model in which an accretion disk is 
broken up by the magnetic fleld of the primary at an in- 
ner radius rj > At rj, matter couples onto the field 
and falls along the field lines towards both magnetic poles 
of the white dwarf. We assume that there is no viscous 
energy release inside rj and that all gravitationally re- 
leased energy is radiated in the flow behind the strong 
shock which forms near the white dwarf surface. This is 
the usual assumption, which allows to treat the guided 
motion as a free fall. It neglects the fact that the ion- 
ized matter in the funnel has to get rid of the particle 
momenta perpendicular to the field and docs so by radia- 
tive losses. Nevertheless, photoabsorption of X-rays in the 
guided stream and reemission of the absorbed energy is 
in line with the model and may cause the funnel to be 
an intense source of UV/optical radiation. Outside n, we 
assume a circularly symmetric steady state disk in which 
the individual annuli radiate at local effective tempera- 
tures Teff(r) as given by Shakura & Sunyaev (1973). In 
addition, we include an inner ring at n, which mimics the 
boundary layer and emits the energy released when the 
matter is braked from the Kepler velocity to the angular 
velocity of the white dwarf. Following Bonnet-Bidaud et 
al. (1982) and Mouchct (1983), we use blackbody spectra, 
which has the advantage of simplicity, but the disadvan- 
tage of systematically distorting the spectrum, especially 
at short wavelengths. For the outer disk, we allow for the 
possibility of heating by radiation from the central source 
or from the elevated magnetically guided flow by letting 
Tcff(r) not drop below a temperature T^i^. This is ob- 
viously a very rough approach, which we justify by its 
simplicity. 

The location of the inner edge of the disk is restricted 
by two requirements. As a first condition, n should not 
exceed the circularization radius rdrc of matter with the 
specific angular momentum it carries over from the inner 
Lagrangian point. No disk can form for n > rcirc- For 
Ml = 0.93 - I.I8M0 and M2 = 0.25 - 0.40 Mg, rare 
falls in the range of 1.06 10^° to 1.48 10^° cm. The sec- 
ond condition results from the observation that the white 
dwarf in V1223Sgr is being spun down, which indicates 
that it should rotate close to equilibrium. The equilibrium 
inner radius of the disk rj c is located inside the corota- 
tion radius r^o = {GMiP^ /An^f/^ = 1.23 lO^^My-'^cm, 
with Ps — 745.5 s the spin period of the white dwarf in 
V1223Sgr (Osborne et al. 1985). The theory of Ghosh 
& Lamb (1979) with later corrections by Wang (1987) 
predicts n^e — 0.98 fco- This does not account, however, 
for the synchronization torque between primary and sec- 
ondary star. Warner (1996) estimates this torque from 
considerations of the magnetism of the secondary star. For 
conciseness, we refer to his Eqs. (22), (23), and (24), and 
note that the magnitude of the effect depends not only on 
the uncertain magnetic moment of the secondary star, but 
also on the ratio between inner edge and spherical Alfven 
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radius, usually taken to be rj/r^ sph — 0.52, and on the 
poorly known screening of the magnetic moment of the 
primary by the accretion disk. In our preferred model B 
(sec below), the disk in V1223Sgr has a rather large cen- 
tral hole and may be quite flimsy. We assume, therefore, 
that screening is rather weak {5 ~ 0.2 in Warner's termi- 
nology). For the parameters of V1223Sgr, we then obtain 
Tie — 0.83 

^CO 

1.0 lO^" cm, which is smaller than r^irc 
for the entire range of possible values of Mi and M2. 

This theory also allows to estimate the mag- 
netic moment of the primary and its surface field 
strength. We use 0.52rju,sph with r^,sph = 

2-3/7(GMi)-i/^M-2/7^4/7(,j^^ ^Yiere ^ = BiRf is the 

magnetic moment, Bi is the surface field strength of the 
white dwarf in the orbital plane, and the dipolar magnetic 
axis is assumed to be aligned with the rotational axis. In 
the aligned case, the polar field strength is Bp = 2Bi, 
while for an inclined magnetic rotator, the factor between 
Bp and Bi is smaller than two. 

3.5. Derivation of Mi, M^, and 

The free fall of matter from ri to Ri releases an accretion 
luminosity Lacc,i in a shock with temperature kTg 

Lacci = M-GMi(l/i?i-l/n) (6) 
kT, = (3/imu/8) • GMi{l/Ri - l/n), (7) 

where /x = 0.617 is the molecular weight of all particles 

for solar composition, is the unit mass, and prompt 
equilibration between electrons and ions in the shock is 
assumed. X-ray emission is the major channel of energy 
release in the post-shock region, but cyclotron emission 
may also contribute despite the negative results of polar- 
ization measurements (Cropper 1986). We represent the 
cyclotron fraction by a parameter a and write 

Le^d = (1 + a) Lx = (1 + a)(l + /?) Lx,i. (8) 

The quantity a, which depends on the field strength and 
particle density in the post-shock region (e.g. Woelk & 
Beuermann 1996), is of the order of 0.1. 

Eq. (7) can be solved for Mi with kT^ = 43li^keV 
(Beardmore et al. 2000), an adopted value for n and a 
mass radius relation of white dwarfs Ri{Mi). We use the 
relation i?i = [1.463 - 0.885 (Mi/M©)] ■ 10^ cm, which we 
found from fitting the radii of CO white dwarf models by 
Wood (1995) with a thick hydrogen envelope, an effec- 
tive temperature of 3 10^ K, and masses between 0.80 and 
1.20 Mq. We use two values of n, which correspond to 
the two models considered below, n = 8.0 10^ cm and 
n = 1.0 10^° cm (see Sects. 3.6.1 and 3.6.2). The white 
dwarf masses for these two cases are Mi=1.18 ± 0.06 
and Mi=0.93±0.12Mq, respectively. A white dwarf mass 
lower than ~ 0.8 Mq can be excluded, as noted already 
by Beardmore et al. (2000). 

We derive the accretion rate Mx needed to gen- 
erate the emitted X-ray spectrum femiiE) by elim- 
inating GMi(l/i?i-l/ri) from Eq. (6) with help of 



Eq. (7) and replacing Lacc,i from Eq. (8). With _Lx.i 
and /? = 1 from Sect. 3.3.1, we then obtain 
Mx = (1-24 ± 0.24) 10^^(1 -I- a) gs'S where the error de- 
pends primarily on the remaining uncertainty in d and 
the value of a. Note that Mx is independent of the white 
dwarf mass. The X-ray luminosity is Lx = 2Z/x,i = 
(2.23 ±0.71) 10^4 ergs-^ 

3.6. Interpretation of the overall UVptical/IR spectral 
energy distribution 

In what follows, we interprete the observed UV/optical 
spectral energy distribution in terms of two simple mod- 
els and list their parameters in Table 4. Model A is that of 
a luminous accretion disk, which has its inner edge close 
to the white dwarf (Bonnet-Bidaud et al. 1982, Mouchet 
1983). Model B has the inner edge at n = 1.0 10^° cm 
and requires a different origin for most of the UV emis- 
sion. We use a distance of 527 pc and a radial velocity 
amplitude of the white dwarf Ki = 56kms~^ (Section 
3.2). We quote the results in Table 4 for the nominal sec- 
ondary mass M2 = 0.40 Mq and comment on the effect of 
choosing a lower value of M2. 

Figure 2 shows the mean orbital and rotational spec- 
tral energy distribution of V1223Sgr in its normal bright 
mode. It is based on the lUE spectrophotometry and 
optical photometry of Bonnet-Bidaud et al. (1982) and 
Mouchet (1983) (solid circles) and is supplemented by 
the (non-simultaneous) 2MASS infrared photometry from 
Table 1 (open circles). The integrated flux for the wave- 
length range A > 1250 A is given in Table 3, line 5, and an 
estimate of the flux between the Lyman edge and 1250A 
in line 4. 

3.6.1. Model A: steady state accretion disk 

The model of a luminous accretion disk (Bonnet-Bidaud 
et al. 1982, Mouchet 1983) can account for the observed 
UV/optical/infrared emission. The model requires a high 
white dwarf mass, an outer disk radius of 90% of the 
Roche radius, an inner radius close to the white dwarf, 
and an accretion rate Mbb-disk = 3.1 lO^'^gs"^, which 
exceeds Mx from above. The best fit to the spectral flux 
distribution is shown in Fig. 2 (upper data set and curve) 
and the parameters are listed in Table 4. The spacing be- 
tween n and i?i is tailored to yield kTg = 43 keV. The de- 
rived accretion rate Mbb-disk represents an overestimate 
because the blackbody disk and its boundary layer at ri 
emit 2/3 of the total flux in the Lyman continuum, while 
typical CV disks have the Lyman jump strongly in ab- 
sorption. The blackbody disk is a coarse approximation 
and we quote Mbb-disk in Table 4 in brackets to indicate 
this uncertainty. If the secondary mass is reduced from 
0.40 to 0.25 Mq, the blackbody accretion rate increases to 
Afbb-disk = 5.0 10^''gs~^, partly because the now larger 
inclination of 44° implies a smaller projected area of the 
disk. 
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Table 4. Parameters of disk models A and B. For both 
models, we assume M2 = 0.40 M©, Ki = 56kms^^, and 
d = 527 pc (see text). 



1000 2000 5000 10000 20000 

Wavelength (A) 

Fig. 2. Overall energy distribution of V1223 Sgr for A > 
1000 A (see text). The data are shown twice with two dif- 
ferent model fits. Model A (top data set, fluxes multiplied 
by 10) is for a luminous steady state accretion disk and 
model B (bottom data set, fluxes as observed) is for a 
truncated disk plus reprocessed component. 

The real problem of the disk model is the location of 
the inner edge of the disk far inside the corotation ra- 
dius, which is incompatible with the observed spin down 
of the white dwarf (van Amerongen et al. 1987, Jablonski 
& Steiner 1987). The theory presented above suggests that 
Tj should be about an order of magnitude larger than as- 
sumed in this model. Whatever the remaining uncertain- 
ties in the theory, the inner edge of the disk can not be 
close to the white dwarf. 



3.6.2. Model B: truncated disk and reverberation of 
X-rays 

The alternative model B assumes a disk, which is trun- 
cated at n = 0.83 rco = 1.0 10^° cm. Since the truncated 
disk is intrinsically faint, we assume that heating by irra- 
diation raises its flux to the observed level, which requires 
Tirr = llOOOK. The truncated disk is represented in Fig. 2 
by an llOOOK blackbody (short dashed curve). In Table 4, 
section (4), lines 1 and 2, we list its UV/optical flux. The 
contribution from the secondary star is minute in com- 
parison (dotted curve for a 3100 K blackbody adjusted to 
K=15.1). Heating the disk to llOOOK requires that it is 
moderately inflated. For a point source at the white dwarf, 
a solid angle of f2 ~ 2 sr or a half opening angle 7 = 9° is 
needed. If the disk is irradiated by the funnel emission, 7 



Quantity 


Model A 


Model B 


(1 J System, pctvameters: 






White dwarf mass Mi (Mq) 


1.18 


0.93 


White dwarf radius Ri (cm) 


4.17 10^ 


6.42 10® 


Inclination i 


27.6° 


24.4° 


Separation a (cm) 


9.21 10"^ 


8.69 10^° 


Corotation radius Tco (cm) 


1.30 10^° 


1.20 10^° 


Inner disk radius n (cm) 


8.00 10® 


1.00 10^° 


Outer disk radius ro (cm), (0.90rRoche,wd) 


3.48 10^° 


3.21 10^° 


(2) Temperatures: 






Un-irradiated disk: Tmax (K) 


58030 


(6534) 


T^in (K) 


6744 


(4556) 


Irradiated disk: Ti^r (K) 




11000 


'Reprocessed' component (K) 




40000 


(3) Accretion rate and accretion luminosity 


for a = 1:^ 




X-ray accretion rate Mx (gs~^) 


1.24 10^^ 


1.24 10" 



Blackbody accretion rate Mbb-disk (gs ) (3.1 lO'^ ) 

Total accretion luminosity Lace (ergs~^) 

X-ray luminosity Lx (ergs"^) 2.23 10^" 

(4-) Disk and non-disk/reprocessed fluxes at the Earth:^ 

Disk > I250A (ergcm-^s-^) 7.61 10"^° 

Disk 912 - I250A (ergcm-^s-^) 2.85 10"^° 

'Reprocessed' > 1250A (crgcm^^s^^) 
'Reprocessed' 912 - 1250A (ergcm^^s^^) 

(5) Disk and non-disk/reprocessed luminosities: 

Disk > I250A (ergs^i) 1.23 10*" 

Disk 912 - I250A (ergs"^) 4.59 10** 
'Reprocessed' > 1250A (ergs"^) 
'Reprocessed' 912 - 1250A (ergs"^) 



2.34 10** 
2.23 10** 



2.70 10"^° 
0.02 10~i° 
4.88 10"^° 
2.96 10"" 



4.15 10''* 
0.03 10** 

- 8.9 10** 

- 5.4 10** 



^ Mx and Lx refer to both poles of the white dwarf. Lace ~ 
GMiM /Ri. ^ Fluxes corrected for interstellar extinction with 
= 0.47. 

could be smaller. In this model, much of the observed UV 
radiation must be due to the reverberation of X-rays. 

That the observed ultraviolet radiation is not entirely 
of disk origin was indicated already by Mouchet's (1983) 
observation that the UV flux varied by 23% over 1.5 h 
while the V magnitude stayed constant. Furthermore, 
Bonnet-Bidaud et al. (1982, see also Welsh & Martell 
1996) showed that the pulsed component has a steeper 
spectrum than the mean light. When corrected for extinc- 
tion, it is almost as steep as a Rayleigh-Jeans spectrum, 
fi, oc A~^'^. The wavelength dependence of this component 
is depicted by the open triangles in Fig. 2, adjusted to 30% 
of the total mean flux in the B-band. This percentage is at 
the upper end of the observed range of amplitudes (Steiner 
et al. 1981, Bonnet-Bidaud et al. 1982, King & Williams 
1983, Warner & Cropper 1984) and is a reasonable choice 
if the observed variations in amplitude represent variations 
in the visibility of the modulation rather than true varia- 
tions of the pulsed flux. For illustrative purposes, we have 
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adjusted a 40000 K blackbody to fit these fluxes (Fig. 2, 
long dashed curve) and the sum of the two components to 
fit the total observed flux. The 'reprocessed' flux is listed 
in Table 4, section (4), lines 3 and 4. We do not contend 
that the 794 s pulsed component extends into the UV as in- 
dicated by the 40000 K blackbody, but suggest that there 
may exist more than one component of reprocessed light. 
Combined they may account for much of the UV flux. In 
view of Penning's (1985) discovery of radial velocity vari- 
ations at the spin period of the white dwarf, we suspect 
that a spin-modulated component may still be hidden in 
the UV. The heated pole cap of the white dwarf is a viable 
contender. The fact that such a component has not been 
found at optical wavelengths may be due to a small ampli- 
tude of modulation and/or a contrived geometry (Hellier, 
2003). The sum of truncated disk and 'reprocessed' com- 
ponent is seen to fit the data as well as the model of a 
luminous disk. The present model B has the advantage of 
being internally consistent with th(^ spin down of the white 
dwarf and with the X-ray derived accretion rate M-^. Its 
parameters depend only weakly on the choice of M2. 

Is the generation of much of the UV flux by reverber- 
ation of X-rays energetically feasible? In order to test this 
hypothesis, we convert the components of the observed 
flux F in Table 4, section (4), to luminosities L = gcPF 
using appropriate (i.e. non-blackbody) geometry factors 
g. For the 'disk' flux we assume a grey limb darkening 
law, oc (1 + 1.5cosi) cosi, which yields g = 1.85 tt for 
i = 24.3°. For the 'reprocessed' component, g is similar 
to that of the disk if the heated polar cap is the main 
source, while an estimated g ~ 2.6 tt applies to the toroidal 
geometry of the magnetically guided flow inside ri. We 
use an average g = 2.2 n. The resulting luminosities are 
listed in Table 4, section (5). After corrections for the con- 
tributions by the un-illuminated disk and the secondary 
star, we obtain an upper limit to the reprocessed lumi- 
nosity Lrep,max(> 9121) = 1.78 lO^^ergs"! = 0.80 Lx- 
The uncertainty in this number is fairly large because of 
possible errors in the geometry factors and the remaining 
uncertainty in the X-ray luminosity. Long-term variabil- 
ity of V1223 Sgr could affect the result, but is probably 
of minor importance because the AAVSO records (Mattel 
2003) show V1223Sgr at ^ = 13.1 at the time of the 
ASCA X-ray observations, in agreement with its long- 
term mean magnitude. The load on the X-ray source is 
reduced if part of the optical emission is of cyclotron ori- 
gin, represented by the quantity a introduced in Sect. 3.4, 
which depends on the field strength and the pre-shock 
mass flow density rh (in gcm^^s^^). From Mx and an es- 
timate of the area of the ring-shaped accretion region on 
the white dwarf we expect rh 3 — lOgcm^^s^^. For 
B 2i 10 - 15 MG, we then estimate a ~ 0.05 - 0.15 
using the radiation-hydrodynamic calculations of Woelk 
& Beuermann (1996, their Fig. 9, bottom panel, lower 
curve). For a ~ 0.10, L^ep/Lx drops to a comfortable 0.69, 
which can easily be accounted for by X-ray heating of the 
pole cap of the white dwarf and reprocessing in the funnel, 
in the accretion disk, and the irradiated face of the sec- 



ondary star (see Beardmore et al. 2000 for an analysis of 
the energy-dependent X-ray albedo from the white dwarf). 
For a = 0.1 ± 0.05, the accretion rate and the total accre- 
tion luminosity increase over the values quoted in Sect. 3.5 
and in Table 4 to M = (l+a)Mx = (1.36±0.27) lO^^gs'^ 
and Lace = 1.049 (1 + a) Lx = (2.57 ± 0.83) lO^^ergs"!, 
where the errors are from Sect. 3.5 and the numerical fac- 
tor in the last relation accounts for the energy released 
outside Vi. Since the field estimate given below suggests 
that some cyclotron emission should be present, we ac- 
cept the latter values as our best estimates. 

3. 7. Magnetic moment and field strength of the white 
dwarf 

The different sizes of the central holes in the accretion 
disks imply different surface field strengths Bi of the white 
dwarfs in models A and B. A rough estimate of Bi may be 
obtained by equating n with 0.52rp,sph (sec above). For 
model B with n ~ 1.0 10^° cm , /x = 2.06 lO^^ Q^m^, and 
for n = 6.42 10® cm (Table 4), Bi = 8MG. Depending 
on the obliquity of the dipole, the polar field strength is 
in the range of Bp = 10 — 16 MG. The small inner hole 
of the disk with n = 8.0 10® cm in model A implies a 
much smaller surface field strength of about 0.5 MG. Our 
favored model B places VI 223 Sgr at the lower range of 
field strengths observed in polars. Since low-field polars 
are characterized by weak soft X-ray emission, the lack of 
a strong soft X-ray source in V1223 Sgr is not surprising. 

4. Conclusion 

We have presented an accurate parallax of V1223 Sgr 
which allows us to derive the luminosities in the different 
wavelength bands. Based on this result, we have tested 
the hypothesis that much of the UV/optical emission is 
produced by the reverberation of X-rays and not by the 
release of gravitational energy in a luminous disk. Our 
analysis is based on the assumption that we see only the 
X-rays from one pole. This assumption and the high X- 
ray temperature reported by Beardmore et al. (2000) lead 
to a total X-ray luminosity ^ (2.2 ± 0.7) 10^'', which 
is much higher than previously thought and can power 
a large fraction of the observed UV/optical radiation by 
the reprocessing of X-rays. Interestingly, this is possible 
without the presence of a substantial source of XUV ra- 
diation. Likely reprocessing sites are the pole cap of the 
white dwarf, the magnetically guided accretion flow, the 
disk, and the irradiated face of the secondary star. While 
the latter two sites may be responsible for the observed 
flux modulated at the sideband frequency (Steiner et al. 
1981), reprocessing in the former locations may produce a 
so far undiscovered component, which is photometrically 
modulated at the spin period of the white dwarf. The ac- 
cretion rate derived from the X-ray luminosity and X-ray 
temperature is M = (1.4 ±0.3) 10^''gs~^, which includes 
a small correction for the contribution by cyclotron emis- 
sion. This result is independent of the white dwarf mass. 
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The observed UV/optical/IR spectral energy distribu- 
tion can equally well be fitted by a truncated disk plus 
reprocessed component (model B) and by a luminous disk 
(model A). The observed spin-down of the white dwarf 
(Jablonski & Steiner 1987, van Amerongen et al. 1987), 
however, requires that the the inner edge of the disk is 
not too far inside the corotation radius, a condition which 
is met only by the truncated-disk model B. 

Strong internal absorption and the reprocessing of a 
major fraction of the emitted X-rays is a general feature 
observed also in other intermediate polars, some of which 
are even more strongly internally absorbed than V1223 Sgr 
(e.g. Norton & Watson 1989). Simultaneous X-ray/UV 
studies could shed light on the physical processes acting 
in these systems. It is surprising that no such study is yet 
available. 
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